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^ ; Abstract 

We consider integrable models, or in general any model defined by an i?-matrix, 
>• . on random surfaces, which are discretized using random Manhattan lattices. The 

^ i set of random Manhattan lattices is defined as the set dual to the lattice random 

Q\ \ surfaces embedded on a regular d-dimensional lattice. They can also be associ- 

^ ■ ated with the random graphs of multiparticle scattering nodes. As an example 

^ ■ we formulate a random matrix model where the partition function reproduces 

^r^ . the annealed average of the XXZ Heisenberg model over all random Manhattan 

lattices. A technique is presented which reduces the random matrix integration 
in partition function to an integration over their eigenvalues. 

X 



1 Introduction 

One of the major goals of non-critical string theory was to describe the non- 
perturbative physics of non-Abelian gauge fields in two, three and four dimen- 
sions. The asymptotic freedom of these theories allowed us to understand the 
scattering observed at high energies [1, 2], but it also made the long distance, low 
energy sector of the theories non-perturbative and indicated a non-trivial struc- 
ture of the vacuum [3, 4]. It became necessary to develop non-perturbative tools 
which would allow us to study phenomena associated with e.g. confinement. One 
possibility which attracted a lot of attention was the attempt by Polyakov to re- 
formulate the non-Abelian gauge theories as a string theory. This line of research 
led Polyakov to his seminal work on non-critical string theory [5] . He showed that 
the presence of the conformal anomaly forces us to include the conformal factor in 
the string path integral, and that the action associated with the conformal factor 
is the Liouville action. In his approach the study of non-critical string theory 
becomes equivalent to the study of two-dimensional quantum gravity (governed 
by the Liouville theory) coupled to certain conformal matter fields. 

Attempts to understand and define rigorously the quantum Liouville theory 
triggered the lattice formulation, presenting the two-dimensional random sur- 
faces appearing in the string path integral as a sum over triangulated piecewise 
linear surfaces [6, 7, 9]. This sum over "random triangulations" (or "dynamical 
triangulations" (DT)) could be represented by matrix integrals and in this way 
certain matrix integrals became almost synonymous to non-critical string theory. 
Somewhat surprisingly many of the lattice models were exactly solvable and at 
the same time it was possible to solve the continuum quantum Liouville model 
via conformal bootstrap, and whenever results of the two non-pertubative meth- 
ods could be compared, agreement was found. However, the solutions only made 
sense for matter fields with a central charge c < 1 [12]. Thus, in a certain way 
the approach of Polyakov came to a stalemate with regards to the question of us- 
ing non-critical string theory to understand the non-perturbative aspects of three 
and four-dimensional non-Abelian gauge theories, which seemingly require c > 1. 
It is thus of great interest to try to study new classes of random surface models 
which might allow us to penetrate the c = 1 barrier. This is one of the main 
motivations of this paper. We propose to consider a new class of random lattices, 
the so-called random Manhattan lattices. One is led to such lattices by studying 
the random surface representation of the 3d Ising model on a regular 3d lattice 
[13, 14], and via the study of the Chalker-Coddington network model [16]. The 
study of the latter model led to the idea that an i?-matrix could be associated 
to a random Manhattan lattice, and we will consider how to couple in general a 
matter system defined by an i?-matrix to a random lattice. By summing over the 
random lattices (i.e. taking the annealed average) we thus introduce a coupling 




Figure 1: A random surface on the 3d cubic lattice and the construction of the 
dual lattice surface, as described in the text. 

between the integrable model and two-dimensional quantum gravity. 

More precisely we start with an integrable model on a 2d square lattice, as- 
suming we know the R-matrix. We then show that the same R-matrix can be 
used on a so-called Random Manhattan Lattice (RML) (see Fig. 1), which is 
a lattice where the links have fixed arrows which indicate the allowed fermion 
hopping. No hopping is allowed in directions opposite to arrows. The summa- 
tion over the RMLs can be performed by a certain matrix integral related to the 
R-matrix. This matrix integral is somewhat different from the the conventional 
matrix integrals used to describe conformal field theories with c < 1 coupled to 
2d quantum gravity, and thus there is hope than one can penetrate to c = 1 
barrier. Below we describe the construction in detail. 

2 The model 



As mentioned one arrives in a natural way to a RML from the study of the 3d 
Ising model on a regular cubic lattice. The high temperature expansion of the 
Ising model can be expressed as a sum over random lattice surfaces of the kind 
shown in Fig. 1, and on these two-dimensional lattice surfaces one constructs a 
kind of dual lattice by the following procedure: The lattice surface consists of 



plaquettes. Consider the mid-points of the hnks on the plaquettes as sites of the 
dual lattice, and consider arrows on the links as shown in Fig. 2. 





(a) (b) 

Figure 2: Assignment of arrows to dual lattice 

The allover orientation of arrows on the plaquettes should be such that the flow 
to neighbouring plaquettes is continuous as illustrated in Fig. 3. This type of dual 
lattice with arrows will be a finite Manhattan lattice corresponding to a particular 
plaquette lattice surface on the regular three-dimensional lattice. There is a one 
to one correspondence between the plaquette surfaces on the regular lattice and 
the finite Manhattan lattices described above. 

A second way of obtaining a RML is by starting from oriented double line 
graphs, like the ones introduced by 't Hooft, and then modify the double line 
propagator like shown in Fig. 4. 

We will now attach an i?-matrix of an integrable model to the squares of 
the RML with the index assignment shown in Fig. 5. Two neighbouring squares 
will share one of indices, and the same is thus the case for the corresponding 
i?-matrices, and a summation over values of the indices are understood, resulting 
in a matrix-like multiplication of i?-matrices. To a RML Q we now associate the 
partition function 

z{n) = l[R, (1) 

where the summation over indices is dictated by the lattice. Our final partition 
function is defined by summing over all possible (connected) lattices Q: 



J2zm 



-tJ.\n\ 



(2) 



where /i is a "cosmological" constant which monitors the typical size |f2| of the 
lattice Q. As long as we restrict the topology of lattices Q entering in the sum (2), 
there will exist a critical fic such that the sum in (2) is convergent for fi > fic and 
divergent for fi < Hc- We will be interested in a limit where the average value of 
Q becomes infinite, and this limit is obtained by approaching ^c from above. The 
summation over the elements in Q, i.e. the summation over a certain set of random 




Figure 3: Random Manhattan lattice 
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Figure 4: Manhattan lattices and double line graphs 



2d lattices, is a regularized version of the sum over 2d geometries precisely in the 
same way as in ordinary DT. It is known the sum over random polygons (triangles, 
squares, pentagons etc) with positive weights under quite general conditions leads 
to the correct continuum limit, i.e. the functional integral over 2d geometries, 
when the link length goes to zero. Thus it is natural to assume that the sum 
over RML will also represent in correct way the sum over 2d geometries when the 
link length goes to zero. Under this assumption we have coupled a given two- 
dimensional model, integrable on a regular lattice, to two-dimensional quantum 
gravity. In the next Section we will, in order to make the discussion more explicit, 
consider the XXZ Heisenberg model, where the R-matrix is known. 
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Figure 5: Index assignment of the i?- matrix 

3 The matrix model 

In order to represent the XXZ model as a matrix model which at the same time 
will offer us a topological expansion of the surfaces spanned by the oriented ribbon 
graphs considered above, we consider the set of 2N x 2N normal matrices. We 
label the entries of the matrices as Map^ij, where a, (3 takes values 0,1 and i, j takes 
values 1, . . . , A^. The a, /3 indices refer to the XXZ model, while the i, j indices 
will be used to monitor the topological expansion. A normal matrix is a matrix 
with complex entries which can be diagonalized by a unitary transformation, i.e. 
for a given normal matrix M there exists a decomposition 

M = UM^'^^U^ (3) 

where f/ is a unitary 2N x 2N matrix and M^^^ a diagonal matrix with eigenvalues 
mjj ^j which are complex numbers. 
Consider now the action 

S{M) = M:^^^fi-jM,,^,,, - E f t^ (^" + (^')") • (4) 

n=3 

We denote the sum over traces of M and iVf^ as the potential. The matrix 
partition function is defined by 

'"dMe-^^(^). (5) 

When one expands the exponential of the potential terms and carries out the 
remaining Gaussian integral one will generate all graphs of the kind discussed 
above, with the i?-matrices attached to the graphs as described. The only dif- 
ference is that the graphs will be ordered topologically such that the surfaces 
associated with the ribbon graphs appear with a weight N^, where x is the Euler 
characteristics of the surface. If we are only interested in connected surfaces we 
should use as the partition function 

F = log Z. (6) 



In particular the so-called large A^ limit, which selects connected surfaces with 
maximal X) "^iH sum over to the planar (connected) surfaces generated by F , 
since these are the connected surfaces with the largest x- 

Explicitly, in the case of the XXZ Heisenberg model the R matrix is given by: 

^ ^' — ' ' \/3' , f^ W ^ / N/3' 



:^ (^i). ® (^i);; + MS ® M^ . (7) 
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As an abbreviation we will write 

^.r =^a® <yJa (8) 

where a summation over index a is understood, o"o = 1, the identity matrix, and 
(Jq, a = 1, 2, 3 are the Pauli matrices. 

Our aim is to decompose the integration over the matrix entries of M into their 
radial part Md and the angular [/-parameters. This decomposition is standard, 
the Jacobian is the so-called Vandermonde determinant (also in the case of normal 
matrices, [17]) When we make that decomposition the potential will only depend 
on the eigenvalues nv^ -^ and for the measure we have: 

dM = df/ n d^s, dmit n I (^s. - ^s.) f • (9) 

a,i a,i^l3,j 

However, the problem compared to a standard matrix integral is that the matrices 
U, introduced by the transformation (3), will appear quartic in the action (4). 
Thus the L'^-integration does not reduce to an independent factor, decoupled from 
the rest. Neither is it of the Itzykson-Zuber-Charish-Chandra type. 

In order to perform the integral we pass from the transformation (3) which 
is given in the fundamental representation, to a form where we use the adjoint 
representation. Let us choose a basis t^ for Lie algebra of the unitary group 
U{2N) in the fundamental representation. The normal matrix M can also be 
expended in this basis: 

M = CAt'^, trtV = (5^^, (10) 

where the last condition just is a convenient normalization. For a given U be- 
longing to the fundamental representation of U{2N) the corresponding matrix in 
the adjoint representation, A(f/), and the transformation (3) are given by 

AiU)AB = tTt^Ut''U\ CA = AABCil\ (11) 



where C^ denotes the coordinates of M^^'' in the decomposition (10). The trans- 
formation (3) is now hnear in the adjoint matrix A and the action (4) will be 
quadratic in A. However, we pay of course a price, namely that the entries of the 
(2A^)^ X {2NY unitary matrix A satisfy more complicated constraints than those 
satisfied by the entries of the 2N x 2N unitary matrix U . We will deal with the 
this problem below. First we express the action (4) in terms of the eigenvalues 
mi „• and A. 
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Figure 6: Index assignment of the i?- matrix 



For this purpose it is convenient to pass from the representation (7) of the 
i?- matrix to the cross channel (see Fig. 6), which is achieved by: 



/3/3' 









(12) 



which amounts to making a particle-hole transformation -H- 1 for the indices a' 
and /3. After some algebra one obtains: 



/-^N^/3' 6 + c h — c 

^ [1 ® 1 + fTa ® as] • 



(13) 



We can now write 



(^1Z'=^«®^«^^ 



(14) 



where the summation is over a, a = 0, 1, 2, 3, ctq = 1, the identity 2x2 matrix. 

One has: 

'ab + cb — ca\ ,. _, 

(15) 



(16) 



_2' 2 ' 2 '2, 
Using the cross channel i?-matrix (14) the action (4) becomes: 



Let T^, /i = 1, . . . , N"^ denote generators of the Lie algebra of U{N), appro- 
priately normalized such that 



SI 5^ 



'ij'i'j'- 



(17) 



If we insert (17) into the action (16) we obtain 

S = tT (mVV^) rtr (aV^M) (18) 

where we can view (x°'t'^, a = 0, 1, 2, 3 and fi = 1, . . . , N"^ as the {2N)'^ generators 
t^ of the Lie algebra of U{2N). Formulas (10) and (11) with the generators a"-T^ 
read 

,W = itr (M('^VV^) , (19) 



"^a,M 



\afi,a' fi' 



itr (a'^T^Ua'^'T^'U^) . (20) 



We now want to use (3) and (19) and (20) to express the matrix M in the action 
(18) in terms of ma,l and A'^^'"''^' , and we obtain 

S = mifJ^:A"'^'"'^'rA"'''""^"mJ?^„ - V{mi^^) (21) 

It is convenient to choose following basis for generators Ia, A = {a,ij) = 

1 ■ ■ ■4N'^ of U{2N). For Cartan sub-algebra we take to,ii = 1 §>> 1 for common 
phase factor, ts,ii = cr3^T^\ i = 1- ■ ■ N and to,M = 1(1 + 0-3) (8)(r*'* — r*"^'*"^), i = 

2 ■ ■ ■ N for other 2N — 1 traceless generators, where 

{r^^)^^ = 5,,,6,i. (22) 

For the remaining, non-diagonal generators we take o'aT^-' , a = 1,2, i,j = l---N 
and aaT^^ , a = 0, 1, 2, 3, i ^ j = 1- ■ ■ N . Then, for this choice of generators and 
from ^Jjjj = |tr (^M^'^'Ha^ij) we have 

<Ji = trM = 5^(mH,,, + mJf4j (23) 

rr^^ = rr^i,u-rr^2,u, (24) 

^oji = ^ali,ii - rnl:^li,i-u-i (25) 

rr^l,^^ = mi%^u = (26) 

with the rest of elements m^ajj = 0; a = 0, 1, 2, 3, i ^ j = 1,- ■ ■ N. 

3.1 The space of integration 

We now change the integration over the unitary matrices U{2N) in formula (9), 
which are in the fundamental representation, to the unitary matrices A(2A^) in 
the adjoined representation. Rather than using the Haar measure expressed in 
terms of the ^/-matrices we should express the Haar measure in terms of the 
A-matrices. 
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Since normal matrices M can be regarded as elements in the algebra u(2N), 
the action of A, defined by the formula (11) on its diagonalized form (23), will 
form an orbit in the algebra with the basis consisting of all diagonal matrices. 
Diagonalized elements of M are invariant under the action of the maximal abelian 
(Cartan) subgroup ®f/(l)^^of U{2N). Therefore the orbits are isomorphic to the 
factor space 



U{2N) 



(7(1)(8)-®C/(1)- 

Moreover this factor space is isomorphic to a so-called fiag-manifold, defined 
as follows (see [19, 20] and references there): A single fiag is a sequence of nested 
complex subspaces in a complex vector space C„ 

{0} = Co C a, C ■ ■ ■ C Ca, CCn = C^ (27) 

with complex dimensions dimcCi = i. For a fixed set of integers (ai, 02 ■ ■ ■ a^, n) 
the collection of all fiags forms a manifold, which called a fiag manifold F{ni,n2, ■ ■ ■ Uk) 
where Ui = at — Oj-i. The manifold F{1, 1, ■ ■ ■ 1) is called a full fiag manifold, 
others are partial fiag manifolds. The full fiag manifold F{1, 1, ■ ■ ■ 1) is isomorphic 
to the orbits of the action of the adjoined representation of U{2N) on its algebra 

The set of Cj hyperplanes in Cj+i is isomorphic to the set of complex lines 
in Cj+i. In differential geometry this set is denoted by CP* (and also as Grass- 
manians Gr{l,i)) and is called a complex projective space. Hence, the complex 
projective space is a factor space 

CP^ = ^(^ + ^) = ^ (29) 

U{i)^U{l) f/(l) ^ ' 

where 5*^*+^ is a real 2i + 1 dimensional sphere. 

According to description presented above the orbit of the action of the ad- 
joined representation of U{2N) on the set of normal matrices M is a sequence of 
fiber bundles and locally, on suitable open sets, the elements of the fiag manifold 
can be represented as a direct product of projective spaces (the fibers) 

J^^CP--'xCP--x...CP' (30) 

In simple words we have the following representation of the orbit: any diago- 
nalized normal matrix in the adjoined representation has a following form 



(31) 



V '" V- 

4Af-l 4Af-3 4fc-l 4fc-3 



^3,115 "-") "-") '"0,11 
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The action of the adjoined representation A on this M transforms it into the 
elements of ^^tt,,Jn presented in (30) where CP^'^"^ represents image of the part 



m^'^^ n ■ ■ ■ n 



U{2N) 
®(7(l)2'v 



4fc-l 



This imphes that the measure of our integral over normal matrices M can 
be decomposed into the product of measures of the base space (the diagonal 
matrices) and the flag manifold (the fiber) 



N 



2N-1 



VA = l[dm^^JM% n ^[CP' 



N 



fc=l 

^ C4A:-1 



n^-s^-an^h^ 



i=l 



fe=l 



^4 



V 



s 



4k-3 



S^ 



(32) 



However, since the diagonal matrix elements uIq {^ and rrv^ l^ are complex, our 
action is invariant over ^(1)^^ (one U{1) per marked segment in (31) and we 
can extend the integration measure from (32) to 



N 



N 



VA = \{dmZdrr^l:[{v[S^'-^ V[S'' 



(33) 



j=i 



k=l 



In other words, we suggest that the action of A on the segments rrig ^^, 0, ■ ■ • and 



4fc-l 



M) 



M) 



(d) 



K,kk^ 0, ■ ■ ■ in (31) form vectors m^'zl,^, (r = 1 ■ ■ ■ 4A; - 1) and m\^>z'Q^,^, {s 



4fc-3 



1- ■ ■4:k — 3), respectively, where the real coordinates z^f^ and Zq/^ belong to the 
unite spheres S^''"^ and S^''~^. 

In order to write the measure of integration over the spheres S"^^'^ and S^'^'^ 
we embed them into the Euclidean spaces R^^ and i?^*^"^, respectively, and define 



4fc 



4A; 



p^^4fe-l" 



ik-l 



s=l s=l -^ s=l 

Ak~2 Ak-2 „ 4fc-2 

5( Et^j' - 1) n ^4fc = / d\k n ^^0, 



e-M.(i-E^'l'[-o.J') 



(34) 



where we have introduced Gaussian integrations over the real parameters \a^ki a 



0, 3. These integrations reproduce the factors 



2Vi-Et=T'[<,.P 



which arise from 
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the (5- functions in (34) by integrations over the coordinates z|^ and -ZqI"^- ^^ 
have omitted coefficients ^/vr/2 in front of integrals on the right hand side of the 
expressions (34) since they unimportant for the partition function. 

With this definition of the measure the partition function (5) can be written 
as 

4fc-l 4fe-3 



J 1. n o ^ — 1 ^ — 1 



S{m^aik'^a,k,Za,k) 



A;,a=0,3 s=l s=l 

(35) 
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where W{nv^\^ = W^^^a .• (^ajj — "^p jj)\ is Vandermonde determinant and 



A'' s(a,k) s(a,k) 



s{m%K,,za,,)= Yl [i<lirEi<'^i'^'+^l'^(i-EKJ')-n 



T^a,kk) 



a=0,3,A;=l 6=1 s=l 



Here s(3, k) = Ak — 1 and s(0, k) = 4k — 3, while according to (15) 

I = -{ a,b — c,b + c,a- ■ -b + c,a,a,b — c,b + c,a- ■ ■ a,- ■ ■ a,b — a,b + a, a ] .(37) 

2v^ ^ '^ ^ ' ^ ^ ' ^^v-^y 

4Af-l 4Af-3 3 1 

The length of / precisely is AN"^. The 4 elements a/2, (6 — c)/2, (6 + c)/2, a/2 of 
/^ are placed as in a specific sequence in /, as shown in eq. (37). However, the 
partition function is independent of this choice (which is just our choice arbitrary 
choice) after integration. 

As one can see we have in the partition function (36) simple Gaussian integrals 
over z^f^. These can be evaluated and we are left with integrals over m\lf^ and 
Xa,k only. It is convenient to rescale the Lagrange multipliers and introduce 
^a,k = \'m'bii\~~^'^a,k- Then, after performing the Gaussian integrals, we obtain 

N N , „ N 



a=0,3;fe=l k=l \"''3,kk\ \"''0,kk\ "^ a=0,3;A;=l 

N N 

Yl e-^^^k)-\rr^k^Kk H z^AKk). (38) 



a=0,3;fc=l a=0,3;A;=l 

where 



■^3,fc(A, 



3,fcJ 



ia-xi,r-Hb-c-\yHb+c-\i,r^ 



(« - ^o,k) Hb-c- Xlk) ' {b + c- Ag fc) 2 
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(36) 



Let us demonstrate that the Gaussian integration over the adjoint represen- 
tation matrice A in (38) correctly reproduces the partition function (5) when 
interaction is absent, i.e. V{M) = 0. In this case the integral over normal ma- 
trices M in the fundamental representation of U{2N) can easily be performed 
directly and the result is: 



pr4iV2 



[a\b^ - c^)] 



-Af2/2 



Let us first consider the simple case A^ = L In the general setup this cor- 
responds to having the two shortest, length 3 and 1 segments in the sequence 
(31). The Vandermonde determinant cancels the m's in the denominator in the 
expression (38) of the partition function and integration over the m's leads to 

dA3,idAoj 11 11 

A3,i Ao,i (a - AiJ V2 (5 - c - AiJ V2 ^b + c- AiJ V2 (« _ Ag J V2 

^ (40) 



!(62_c2)]l/2 



\a- 



provided we place the A-poles at zero and the A-branch cuts at different sides of 
the real A-axis. 

Now consider a general N. For a generic i segment in (31) we first represent 
the multipliers Imiii — rnijA and l^^iij — '"^2 jil i^ ^^e Vandermonde determinant 
as a sum over m)j '■■, a = 0, 3: 

i-l 
I (d) _ (rf) I _ I (d) I V^ (d) I 

r=j+l 

I ('^) _ (d-) 1 — 1 W I \^ ( W _ W \\ (A^\ 

l^a=l,M ^a=2,jrjl ~ l^3,M + 2^ V^a=2,rr ^Q=2,r-lr-lJ I • \^^) 

r=j+l 



When using this decomposition in the Vandermonde determinant product, we 
observe that only the contribution of the selected terms Imiil and Imi iJ in 
(41) will cancel all m'^'^^'s in the denominator of (38) and thus lead to a nonzero 
contribution. By Cauchy integration as above we obtain 

N 



[ n i^i^z.,.(A.,.) 

J n Q.._i M.i Ao.j 

(42) 



1 



[a2(?,2 _ c2)]ivv: 
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Other terms in the decomposition (41) will result in Aq^j, a = 0, 3 appearing in the 
denominator of the integral (42) with other powers than one, and the integration 
will give zero for these terms. 

4 Conclusions 

We have defined a matrix model which reproduces the partition function of an 
integrable model (the XXZ model) on random surfaces. The random surfaces 
under consideration appear as random Manhattan lattices, which are dual to 
random surfaces embedded in a rf dimensional regular Euclidean lattice. This 
formulation allows us to consider a new type of non-critical strings with c > 1. 

We have shown that the matrix integral can be reduced to integrals over the 
eigenvalues of matrices. The important ingredient in the integration over angular 
parameters of the matrices (which are usually defined by the Itzykson-Zuber 
integral) is the reduction of the problem to an integration over unitary matrices 
in the adjoined representation. In principle one can now try to apply standard 
large N saddle point methods for solving the resulting integrals over eigenvalues. 

Thoughout in this paper we have considered the Heisenberg XXZ model on 
random surfaces, but the approach can be applied to other integrable models. 
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